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Strangeness-driven phase transition in star matter 
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The phase diagram of a system constituted of neutrons, protons, A-hyperons and electrons is 
evaluated in the mean-field approximation in the complete three-dimensional space given by the 
baryon, lepton and strange charge. It is shown that the phase diagram at sub-saturation densi- 
ties is strongly affected by the electromagnetic interaction, while it is almost independent of the 
electric charge at supra-saturation density. As a consequence, stellar matter under the condition of 
strangeness equilibrium is expected to experience a first as well as a second-order strangeness-driven 
phase transition at high density, while the liquid-gas phase transition is expected to be quenched. 
An RPA calculation indicates that the presence of this critical point might have sizable implications 
for the neutrino propagation in core-collapse supernovae. 

PACS numbers: 26.50.+X, 26.60.-c 21.65.Mn, 64.10.+h, 64.60.Bd, 
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Supernova explosions following the gravitational col- 
lapse of a massive star (M > 8M Q ) are among the most 
fascinating events in the universe as they radiate as much 
energy as the sun is expected to emit over its whole life 
span [l| . Nuclear physics is an essential ingredient in the 
numerical simulations which aim to describe these events, 
since realistic astrophysical descriptions of the collapse 
and post-bounce evolution rely on the accuracy of the im- 
plementation of weak processes and equation(s) of state 
(EOS) 0, H[. Determining the EOS constitutes a par- 
ticularly difficult task since phenomenology ranges from 
a quasi-ideal un-homogeneous gas to strongly interact- 
ing uniform matter and, potentially, deconfined quark- 
matter. The situation is even more difficult if phase 
transitions are experienced, since mean-field models fail 
in such situations 

The Coulomb-quenched liquid-gas (LG) phase transi- 
tion taking place at densities smaller than the nuclear sat- 
uration density (no = 0.16 fm~ 3 ) is, probably, the most 
notorious and best understood case p-Q • At highest den- 
sities, a quark-gluon plasma is expected, but predictions 
on the exact location of the transition are strongly model 
dependent (Io| . In the intermediate density domain sim- 
ple energetic considerations show that additional degrees 
of freedom may be available, such as hyperons, nuclear 
resonances, mesons or muons (TT|. The possibility that 
the onset of hyperons could pass via a first order phase 
transition in neutron stars has been evoked in Ref. 121. 
using a relativistic mean field model, and in Ref. 13j . 
a phase transition between phases with different hyper- 
onic species has been observed for cold matter. Using 
a simple two-component (n, A) model, we have recently 
studied the complete phase diagram of strange baryonic 
matter showing that it exhibits a complex structure with 
first and second order phase transitions (lij . However, 
the exploratory calculation of Ref. [3] neglects the fact 
that in addition to baryon number B and strangeness 



S, the charge Q and lepton L quantum numbers are 
also populated. The thermodynamics of the complete 
system should thus be studied in the four-dimensional 
space of the associated charges Tig , n<j, nz, , n<g . The strict 
electroneutrality constraint uq = 0, necessary to obtain 
a thermodynamic limit 15], makes the physical space 
three-dimensional. As it is known from the EOS stud- 
ies at sub-saturation density [T(|, the introduction of 
the charge degree of freedom can have a very strong in- 
fluence on the phase diagram and cannot be neglected. 
In this work we therefore introduce a four-component 
model constituted of neutrons, protons, electrons and A- 
hyperons. 

In the widely used mean-field approximation (Til ll~7r - 
[23[ the total baryonic energy density is given by the 
sum of the mass, kinetic and potential energy den- 
sity functionals which represents a surface in the three- 
dimensional space defined by the baryon, strange and 
charge density given, in our case, by ub = n n + n p + n\, 
ns = —riA and uq = n p . In the non-relativistic for- 
malism valid in the considered domains of density and 
temperature it reads 



es = 



{ni-rriiC 2 4- ^^ T ^j + e P ot(n n , n p , n A ) (1) 



The single-particle densities are given by the Fermi in- 
tegrals 
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is the Fermi-Dirac inte- 
is the effec- 



whereF,( 7? ) - jQ <~ 1+axp(B _ l?) 
gral, /3 = T^ 1 is the inverse temperature, m 
tive «-particle mass and jli is the effective chemical po- 
tential of the i-species self-defined by the single-particle 
density. 
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A full thermodynamics characterization of the system 
is provided by the pressure P = Tsb — &b + SiA t i n i 
together with the entropy density sb in mean-field, 
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(3) 



The thermodynamical definition rij = f^'J 1/9 allows 
to infer the relation among the chemical potentials /i, 
and the effective parameters jli as /i; = /i, + rn^c 2 + f/,-, 
with J7j = de pot /drii. 

Within the numerical applications we shall use the po- 
tential energy density proposed by Balberg and Gal |24| . 

e P ot(n n ,n p ,n A ) = \n,,ii,n, ■ '/,,/,',//,//, (4) 

ij = {",P,A} 



accounting for nucleon-nucleon, nucleon-A and A-A in- 
teractions, tj denotes the third isospin component of 
particle i. In the non-strange sector the form of the inter- 
action is the same as in the widely used Lattimer-Swesty 
[2o| EOS. The corresponding values of the coupling con- 
stants [24| are reported in Table HI 

Matter stability with respect to phase separation can 
be checked in any point of the extensive variable space 
by analyzing the eigen- values of the curvature matrix 
0, [H S3, Qj = & t f({n i }i={ij,k})/dnidn j , where 
k — B,S,Q and / = e to t — Tstot is the total free- 
energy. The occurrence of, at least, one negative eigen- 
value in a certain domain of (ng, ng, tiq) means that the 
system is unstable versus phase separation. The associ- 
ated 3-dimensional Gibbs construction can be reduced 
to a simpler 1-dimensional Maxwell construction [28[ 
by performing a Legendre transformation with respect 
to two out of the three chemical potentials (1b = fJ-m 
/is = — MA + Mn an d (J<q = n P — (J-n- We have chosen to 
work in the hybrid ensemble (ris, /is, iiq) defined by: 

fbaryon(nB, fJ-S, A*q) = fbaryon ~ MS^S ~ ^Q^Q, (5) 

The upper part of Fig. [1] illustrates the projection of 
the T — phase diagram in the riB — rig plane for /is = 0. 
Two phase-coexistence domains may be identified. The 
one lying along n$ — at sub-saturation density corre- 
sponds to the well known LG like phase transition taking 
place in dilute nuclear matter }28| . The second domain 
lies at supra-saturation densities (tib ^ 2no), and the di- 
rection of phase separation is dominated by the strange 
density. This transition is consistent with our previous 
findings within a simpler 2D model [lH • It is at first sight 
surprising to observe that the coexistence border of the 
latter transition is given by simple straight lines. In prin- 
ciple the coexistence borders of a first order phase tran- 
sition with three conserved charges are given by two sur- 
faces in the three dimensional space. For their projection 
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FIG. 1: (Color online) Borders of the phase-coexistence 
domains at T=0 and us = 0. Upper (middle): 
(n,p, A)-mixture in tib — ns (tib — nc) coordinates. Lower: 
(n,p, A, e)-mixture in ns — ul coordinates. Red: liquid-gas 
phase transition of non-strange dilute nuclear matter; blue: 
non-strange to strange phase transition. The arrows mark 
the directions of phase separation. 



on the tlb — ns plane to be given by a one-dimensional 
curve, these surfaces have to be perpendicular to that 
plane, that is independent of uq. The observed indepen- 
dence on the electric charge shows that the strange charge 
is the dominant order parameter for this transition. How- 
ever, we can expect that some dependence on the electric 
charge would arise if charged strange particles were in- 
cluded, because of the correlation which would then exist 
between uq and ns- The middle part of Fig. [T] illus- 
trates the projection of the phase coexistence domains 
in the ns — nQ plane for /is = 0, corresponding to the 
strangeness-equilibrium condition which is relevant for 
star matter. As discussed above, the coexistence domain 
being three-dimensional this representation depends on 
the value of the third variable given by (is (or ns)- The 
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TABLE I: Coupling constants corresponding to the stiffest interaction proposed in Ref. [2 
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FIG. 2: (Color online) Left: Borders of the strangeness 
driven phase transition domain corresponding to the 
neutral net-charge (n,p, A, e)-mixture at T=10, 20, 30 
MeV and /is = in tib — coordinates. The 
dot-dashed line marks, for T=20 MeV, a path of 
constant Y e — 0.298. Right: Electron fraction, Y e , and 
Hb at the corresponding critical temperature for /is = 0. 



well-known isospin dependence of the LG phase transi- 
tion occurring at n$ = [28[ is apparent. We have just 
noticed that the order parameter of the strangeness phase 
transition is given by a combination of the strange and 
baryonic density. Not surprisingly, the direction of phase 
separation of this transition in the ns — tiq plane is thus 
dominated by the baryonic density. The order parameter 
component in the direction of the electric charge can be 
understood as due to the correlation between the differ- 
ent densities. We are facing a transition between a rel- 
atively diluted, non-strange phase to a relatively dense, 
more strange one. Since A's are neutral, the positively 
charged component of the baryonic density is relatively 
less important in the dense phase, which explains the 
slope of the separation direction. 

We now turn to investigate the influence of Coulomb 
effects on the phase diagram. For simplicity, we will con- 
sider only electrons and neglect other charged leptons or 
mesons. 

Electrons are coupled to charged baryons through 
the electromagnetic interaction, which can modify the 
baryonic phase diagram. However, the charge neutral- 
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ity condition hq = makes the associated chemical 
potential hq ill-defined and keeps the problem three- 
dimensional [15]. Since in homogeneous matter with the 
6 condition ?iq = the Coulomb interaction exactly van- 
ishes [l5T ] , the total mean field pressure can be written as 
a sum over independent terms Pbaryon + Piepton + P 7 + .... 
We shall still concentrate on the Pbaryon contribution, as 
the other terms do not affect the convexity properties 
of the thcrmodynamical potential on which phase tran- 
sitions rely. We have constructed the full phase diagram 
0.2 of the (npAe) system in the (jiB,ns,nT,) space using the 
hybrid ensemble 

fbaryon(nB, fJ-S, = fbaryon — ^S n S — I^L n L ■ (6) 

o In practice, the charge neutrality condition gives n p = n e , 
which allows to infer the electron chemical potential, /i e , 
and thus [II = A*p — Mn + Me- Note that n e = n e - —n e + = 
til(= n p) stands here for the net electron density. 

The lowest panel of Fig. [1] depicts the phase coex- 
istence regions of the (npAe) system at T = and 
/is = 0. In agreement with the results of Ref. |29j. 
a strong Coulomb quenching of the LG-phase transi- 
tion is obtained. However, the coexistence domain of 
the strangeness-driven phase transition is practically un- 
modified. This can be easily understood from the fact 
that the effect of the neutrality condition tiq — on 
the two phase transitions is very different. The phase 
transition at sub-saturation density has the total bary- 
onic density as order parameter. At such densities, ns 
is strongly correlated to n p because of the nuclear sym- 
metry energy which favorizes symmetric n n — n p mat- 
ter. The phase transition thus implies a discontinuity in 
n p = n e , which is strongly disfavored by the huge elec- 
tron incompressibility [16J. At supersaturation densities 
the order parameter is given by ns which is very loosely 
correlated to uq . The phase transition thus does not im- 
ply any strong change in the electron distribution and 
the presence of electrons thus does not influence much 
the phase diagram. 

The temperature dependence of the phase diagram 
along ps = is presented in the left panel of Fig. [5] 
We can observe that the direction of phase separation 
is almost independent of T. More interesting, starting 
from a finite value of T, a critical point appears and sur- 
vives up to very high temperature. On the right panel of 
the same figure the critical temperature and the electron 
fraction Y e = n e /riB are shown as a function of baryon 
density. These values are typically reached within the 
cooling proto-neutron star, meaning that effects of criti- 
cality should be experienced. 
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To explore this issue more quantitatively, we turn to 
calculate the mean free path for the neutrino scattering 
off n, p, and A particles including the long-range corre- 
lations, essential for the study of criticality, in the linear 
response approximation. In the non-relativistic limit for 
the baryonic components, the mean free path at temper- 
ature T of a neutrino with initial energy E v is given by 
1/A = l/\ v + 1/\ A [U [U where the contribution of 
the vector channel is defined as, 



f 



G 



\ V (E„,T) 16 
and that of the axial channel 



^ J(l + cos8)S v (q,T)(l-Mk 3 ))dk 3 , 

(7) 



1 



X A {E^T) 16tt 2 



(3-cos0)5 A (g,r)(l-^(k 3 ))dk3. 

(8) 

In Eqs. ||7J)-||HJ), Gf is the Fermi constant, 8 is the angle 
between the initial and final neutrino momentum (=k3), 
q is the transferred energy- momentum, q = (u>, q), and 
/„ is the Fermi-Dirac distribution of the outgoing neu- 
trino. S v (S A ) are the dynamical response function in 
the vector (axial) channel. Since this study is focused on 
the impact of the density fluctuations close to the critical 
point, only the vector channel is considered. For densi- 
ties close to the critical point, spin-density fluctuations 
are however expected to be small [H, HH . 

The dynamical response function in the vector channel 
is defined as 



S v (q,T) 



1 



7T 1 — CXp(— Uj/T) 



-v 



4 n"(g,r) 



. (9) 



where H v (q,T) is the vector-polarization matrix for the 
three species n, p, and A, given by the Lindhard func- 
tions in the case of the mean-field approximation and by 
the solution of the Bethe-Salpeter equations in the case 
of the mean-field+RPA approximation [3(1 HH, [H| . The 
vector coupling constants are set to be: -1 (n), 0.08 (p), 
-1 (A) 34j]. The residual p-h interaction is derived from 
the potential energy (0) and is closely related to the cur- 
vature matrix without electrons 27]. 

The neutrino mean free path along an arbitrary Y e = 
0.2981 trajectory in the phase diagram which passes by 
the critical point (see dot-dashed line in Fig. [2ja)) is 
shown in Fig. [3] As expected, the RPA correlations 
strongly reduce the neutrino mean free path close to the 
critical point, similar to the critical opalescence effect 
observed for the photon scattering off matter in critical 
water. The ratio of the neutrino mean free path in mean- 
field+RPA approximation over that at the mean-field 
level is shown in panel (b) exploring different neutrino 
energies around the neutrino chemical potential defined 
at beta equilibrium. The effect of the RPA correlations 




FIG. 3: (Color online) (a) Neutrino mean free path for 
the scattering off n, p, and A at T — 20 MeV along a 
constant- Y e = 0.2981 trajectory in the phase diagram 
for E v — \x v , /i^ ± T as a function of the baryonic 
density p. The result of the mean-field approximation is 
compared to the mean-field+RPA. (b) The ratio of the 
mean free path within mean-field+RPA over mean-field 
approximation is shown. 



around the critical point is almost independent of the 
neutrino energy in agreement with the interpretation as 
critical opalescence. 

To conclude, we have studied the phase diagram of 
a mixture constituted of interacting neutrons, protons 
and A-hyperons under the condition of strangeness- 
equilibrium, relevant for supernovae and neutron star 
physics. At supra-saturation densities, a strangeness- 
driven phase transition can take place, depending on 
the assumed strengths of nucleon-A and A-A interactions 
[lij ]. This second transition survives the screening ef- 
fect of electrons and persists over a large domain of tem- 
peratures such that it may have an impact on star phe- 
nomenology. For a first study of this equation of state 
(EoS) within core-collapse supernovae, see [3^. In ad- 
dition to the EoS, linear response theory shows that the 
neutrino mean-free path dramatically decreases close to 
the critical point of this phase transition, which occurs in 
a thermodynamic domain accessible to newly-born proto- 
neutron stars. 

These results present a first step, and quantitative 
results might be somewhat modified in the presence 
of other strange- and non-strange baryonic, leptonic or 
mesonic degrees of freedom. This work is in progress and 
it will make the subject of a forthcoming publication. 
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